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Abstract 

We survey various Ua(1) problems and attempt to resolve the two puzzles 
related to the eta mesons that have experimental verification. Specifically, we 
first explore the Goldstone structure of the r\ and rj mesons in the context of 
rj-rj' mixing using ideas based on QCD. Then we study the eta decays rj — ► 3ir°, 
rj — ► 37T° and rj — > rjmr. Finally we arrive at essentially the same picture in 
the dynamical scheme based on consistently coupled Schwinger-Dyson and 
Bethe-Salpeter integral equations. This chirally well-behaved bound-state 
approach clarifies the distinction between the usual axial-current decay con- 
stants and the 77 decay constants in the 77-7/ complex. Allowing for the effects 
of the SU(3) flavor symmetry breaking in the quark-antiquark annihilation, 
leads to the improved 77-7/ mass matrix. 
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I. INTRODUCTION 



Various statements of "U^(l) problems" related to the eta mesons and their possible 
resolutions have appeared in the literature now for almost three decades. Considerations of 
the eta Ua(1) vacuum Ward identity were discussed by Glashow [l|], Weinberg 0, Crewther 
|| and collaborators. The 11^(1) axial current and its anomalous addition were studied 
by Kogut and Susskind ||]. Semiclassical instantons with topological winding number were 
used by 't Hooft 0. Lastly, the large N c limit together with the 9 vacuum were explored 
by Witten ||. All of the above notions were invoked to resolve the Ua(1) problem. These 
above Ua(1) problems have recently been rekindled by 't Hooft in his text J7J. 

We prefer to focus on two U J 4(l)-type problems that have empirical resolutions and which 
also have a theoretical basis: 

1. Goldstone boson structure of the observed || 7/(547) and 7/'(958) mesons via 77-77' 
mixing in the context of QCD. 

2. Observed eta hadronic decay rates: 

(a) T(r] — > 3n°) = 380 ±36 eV || appears large since it should vanish by the Suther- 
land theorem 0, or be a factor of two smaller in the context of chiral perturbation 
theory |IJ. 

(b) r(j/' 3tc°) = 313 ± 58 eV M appears relatively suppressed because rj' — > 3n° 
phase space is six times larger than for 7/ — > 3tt°. 

(c) F(j]' — > 7/7T7r) = 131 ± 8 keV || is a strong decay, whereas the smaller 3ir decays 
in 2a, 2b above change isospin by one unit and are non-strong decays proceeding 
through the quark mass difference rrid — m u . 

(d) We invoke the AI = 1 w 3 = gA 3 g Coleman-Glashow (CG) |TT| quark tadpole to 
support the current- current Sutherland |§ suppression of the 77 — > 3n decay rates. 
The CG tadpole also explains all 13 hadron (P, V, B, D) SU(2) mass splittings 
|TT| , [12| . Then we use PCAC Consistency |L3| to compute the 7/, 77' — > 3n° decay 



rates in 2a, 2b above. 

The above problems are analyzed in Sees. II and III primarily on the basis of the input 
from meson phenomenology. However, the underlying notions of the quark model are also 
crucial in this analysis. Therefore, in Sec. IV we show the consistency of some of the results 
of Sees. II and III with a sophisticated quark model which has strong and clear connections 
with the fundamental theory - QCD. It is based on the so-called coupled Schwinger- Dyson 
(SD) and Bethe-Salpeter (BS) approach in which one, by solving the SD equation for dressed 
quark propagators of various flavors, explicitly constructs constituent quarks. They in turn 
build qq meson bound states which are solutions of the BS equation employing the dressed 
quark propagator obtained as the solution of the SD equation. If the SD and BS equations 
are so coupled in a consistent approximation, the light pseudoscalar mesons are simultane- 
ously the qq bound states and the (quasi) Goldstone bosons of dynamical chiral symmetry 
breaking (D%SB). The resulting relativistically covariant constituent quark model (such as 
the variant of Ref. ||14|| ) is consistent with current algebra because it incorporates the correct 
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chiral symmetry behavior thanks to D%SB obtained in an, essentially, Nambu-Jona-Lasinio 
fashion, but the former model interaction is less schematic. In Refs. [ 14 — P~9f1 for example, 
it is combined nonperturbative and perturbative gluon exchange; the effective propagator 
function is the sum of the known perturbative QCD contribution and the modeled nonper- 
turbative component. For details, we refer to Refs. ||l4|-|l~8|j , while here we just note that 
the momentum-dependent dynamically generated quark mass functions Aif(q 2 ) (i.e., the 
quark propagator SD solutions for quark flavors /) illustrate well how the coupled SD-BS 
approach provides a modern constituent model which is consistent with perturbative and 
nonperturbative QCD. For example, the perturbative QCD part of the gluon propagator 
leads to the deep Euclidean behaviors of quark propagators (for all flavors) consistent with 
the asymptotic freedom of QCD |L7]. However, what is important in the present paper, is the 
behavior of the same mass functions M. f{q 2 ) for low momenta [q 2 = to — q 2 ~ (400 MeV) 2 ], 
where A4f(q 2 ) (due to D^SB) have values consistent with typical values of the constituent 
mass parameter in constituent quark models. For the (isosymmetric) u- and <i-quarks, our 
concrete model choice |14| gives us .M Uj d(0) = 356 MeV in the chiral limit (i.e., with van- 
ishing m u> d, the explicit chiral symmetry breaking bare mass term in the quark propagator 
SD equation, resulting in vanishing pion mass eigenvalue, = 0, in the BS equation), and 
•Mu,d(0) — 375 MeV [just 5% above M. u ,d(0) i n the chiral limit] with the explicit chiral sym- 
metry breaking bare mass m u ^ = 3.1 MeV, leading to a realistically light pion, = 140.4 
MeV. Similarly, for the s quark, Ai s (0) = 610 MeV. The simple-minded constituent mass 
parameters, denoted below by rh in the case of the isosymmetric u and d quarks, and by m s 
in the case of the s quarks, have therefore close analogues in the coupled SD-BS approach 
which explicitly incorporates some crucial features of QCD, notably D^SB. 



II. GOLDSTONE STRUCTURE OF ETA MESONS 

To resolve Ua(1) problem one, we invoke the U(3) pseudoscalar nonet structure (n, K, rj, 
t]') along with the Gell-Mann-Okubo mass formula m 2 + 3m 2 g = 4m| requiring an octet 
eta mass m m « 567 MeV. While this r] 8 mass is presumed to vanish in the SU(3) x SU(3) 
chiral limit (CL), the companion singlet t)q mass is not expected to vanish in the CL. Using 
the standard relation mixing t], rj' to rj s , rj away from the CL one knows 

m m + m ?o = m \ + m %> ~ L22 Gey2 > or m m ~ 947 MeV (!) 

for masses ^(547), 7/(958), 77 8 (567). 

Here we assumed the standard, most traditional representation of the physical isoscalar 
pseudoscalars rj and 7/ as the orthogonal mixture 

|r/) = cos6» \t] 8 } - sin6> \r} ) , (2a) 
\rf) = sin 9 \rj$) + cos 9 \rj ) , (2b) 

of the respective octet and singlet isospin zero states, r] 8 and r7 . In the flavor SU(3) quark 
model, they are defined through quark-ant iquark (qq) basis states \f f) (/ = u, d, s) as 

|?7 8 ) = -^=(\uu) + \dd) — 2\ss)) , (3a) 
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\ Vo ) = —(\uu) + \dd) + \ss)) . (3b) 

The exact SU(3) flavor symmetry (m u = = m s ) is nevertheless badly broken. It is 
an excellent approximation to assume the exact isospin symmetry (m u = ma), and a good 
approximation to take even the chiral symmetry limit (m„ = = 0) for u and <i-quark, but 
for a realistic description, the strange quark mass m s must be significantly heavier than m u 
and rricL- [In particular, the CL is obviously phenomenologically unrealistic in the strange 
sector, although it is qualitatively meaningful, and in fact useful as a theoretical limit in the 
discussions of the Ua(1) problem.] Thus, with \uu) and \dd) being practically chiral states 
as opposed to a significantly heavier \ss), Eqs. ([|) do not define the octet and singlet states 
of the exact SU(3) flavor symmetry, but the effective octet and singlet states. Hence, as in 
Ref. [0 for example, only in the sense that the same qq states \ff) (/ = u,d,s) appear 
in both Eq. ([2a]) and Eq. ( f2"T5|) do these equations implicitly assume nonet symmetry (as 
pointed out by Gilman and Kauffman ||20|| , following Chanowitz, their Ref. [8]). However, 
in order to avoid the Ua(1) problem, this symmetry must ultimately be broken at least at 
the level of the masses. In particular, it must be broken in such a way that 77 — ^ 773 becomes 
massless but rj' — > 770 remains massive (as in Ref. JT6|) when CL is taken for all three flavors, 
m u ,m d ,m s — > 0. 

Alternatively, one can work in a nonstrange (iV5*)-strange (S) basis \t)ns) an d \vs), where 



\Vns) = + \ dd )) = ^ + y 3 \Vo) , (4a) 

\Vs) = \ss) = -y^|l%) + ^l^o) • (4b) 
In analogy with Eq. (P, in this basis one finds 

m ms + m vs = m l + m % ~ 122 GeV2 ' ( 5 ) 
since (t]'\t]) = and since the NS-S mixing relations are 

\rj) = cos(f)\r] NS ) - sin0 |^ s ) , (6a) 
\rf) = sin 4> \r] NS ) + cos <p \q s ) ■ (6b) 

The more familiar singlet-octet state mixing angle 6, defined by Eqs. ([|), is geometrically 
related to the NS-S state mixing angle (ft above as 

6 = <p- arctan V2 = - 54.74° . (7) 

Although mathematically equivalent to the r]s-T]o basis, the NS-S mixing basis is more 
suitable for some quark model considerations, being more natural in practice when the 
symmetry between the NS and S sectors is broken as described in the preceding passage. 
There is also another important reason to keep in mind the \f]Ns)-\vs) state mixing angle 
0. This is because it offers the quickest way to show the consistency of our procedures 
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and the corresponding results obtained using just one (8 or 0) state mixing angle, with 
the two-mixing-angle scheme considered in Refs. [22-^8|, which is defined with respect to 



the mixing of the decay constants. Namely, in the NS-S basis, the two decay-constant- 
mixing angles in the two-angle scheme (<fi q and S ) are very close to each other. This basis 
thus has the advantage that even in the case of the decay-constant-mixing, the reduction 
to just one mixing angle occurs in a good approximation and one can use just one mixing 
angle p4|- [2"6j , [2lj[ ] around <p q « S . Moreover, phenomenology seems to justify the central 



assumption of Feldmann, Kroll and Stech (FKS) |26| that in the NS-S basis, the decay 
constants follow (in a good approximation) the pattern of particle state mixing, so that the 
NS-S state mixing angle ~ <p q ~ (f) s . On the other hand, when we express our results 
in terms of through the relation (f?]), which always holds since we use a state mixing 
angle, we find they are consistent with the FKS scheme. The application of the two-mixing- 
angle scheme is relegated to the Appendix, where we give our predictions for the 77-7/ decay 
constants, as well as 9 8 and 6*0, the two decay-constant-mixing angles (in the octet-singlet 
basis), in that scheme. Here, we simply note that our considerations will ultimately lead us 
to ~ 42°, practically the same as the result of FKS scheme and theory (e.g., in Table 2 
of Feldmann's review p8|), and in agreement with data. World 77-77' mixing angle data in 
1989 led to [gg] 



41° ±2° or 9 = -U° ±2°. (8) 



A more recent detailed analysis []30| based on 1996 data for decays tensor to pseudoscalar 
T — > PP, radiative vector to pseudoscalar (or vice versa) V — > Pj, P — > Vj, double 
radiative decays 77 — > 77, 77' — > 77, and J/ip — > VP decays (14 such decays) leads to the 
empirical 77-77' mixing angles 

= 43° ±5° or = 42° ±3° (9) 

found respectively from observed branching ratios, B{a 2 — > r]n/KK) = 2.96 ±0.53, B{a,2 — > 
ti'-k/titt) = 0.037 ± 0.007, in complete agreement with (|8|). The 77-77' mixing angles in (H) or 
(0) (for 4 of 14 determinations) depend on the constituent quark mass ratio m s /m ~ 1.45, as 



already found from baryon magnetic moments pl| , meson charge radii |E| and K* — > K^y 
decays El. (771 denotes the isosymmetric average mass m ud .) 

As for a theoretical determination of the 77-77' mixing angle or 9 = — 54.74°, we 
follow the path of Refs. [J21]J . The contribution of the gluon axial anomaly to the singlet 770 
mass is essentially just parameterized and not really calculated, but some useful information 
can be obtained from the isoscalar qq annihilation graphs of which the "diamond" one in 
Fig. 1 is just the simplest example. That is, we can take Fig. 1 in the nonperturbative 
sense, where the two-gluon intermediate "states" represent any even number of gluons when 
forming a C + pseudoscalar qq meson JJT], and where quarks, gluons and vertices can be 
dressed nonperturbatively, and possibly include gluon configurations such as instantons. 
Factorization of the quark propagators in Fig. 1 characterized by the ratio X w m/m s leads 



to the pseudoscalar mass matrix in the NS-S basis |21 



/ ml + 2(3 V2(3X \ / m% \ 

\^[2(3X 2m 2 K -ml + f3X 2 U I ' {LU) 
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where (3 denotes the total annihilation strength of the pseudoscalar qq for the light flavors 
/ = u,d, whereas it is assumed attenuated by a factor X when a ss pseudoscalar appears. 
(The mass matrix in the 773-770 basis reveals that in the X — > 1 limit, the CL-nonvanishing 
singlet t] mass is given by 3(3.) The two parameters on the left-hand-side (LHS) of (|10D , (3 
and X, are determined by the two diagonalized rj and rj' masses on the RHS of (|T0D. The 
trace and determinant of the matrices in ( fLOj) then fix (3 and X to be |21] 



4(m 



0.28 GeV 2 



X « 0.78 



mi 



(11) 



with the latter value suggesting a constituent quark mass ratio X -1 ~ m s /m ~ 1.3 , near 
the values in Refs. |2"9"|-|3~5[] , m s /m m 1.45. 

This fitted nonperturbative scale of (3 in ([□]) depends only on the gross features of QCD. 
If instead one treats the QCD graph of Fig. 1 in the perturbative sense of literally two gluons 
exchanged, then one obtains |3^] only {3 2 g ~ 0.09 GeV 2 , which is about 1/3 of the needed 
scale of (3 found in (|TTD . (This indicates that just the perturbative "diamond" graph can 
hardly represent even the roughest approximation to the effect of the gluon axial anomaly 
operator t f^G^G 1 * v .) The above fitted quark annihilation (nonperturbative) scale f3 in 
([□]) can be converted to the NS-S rj-r]' mixing angle <fi in (H) from the alternative mixing 
relation tan 20 = 2\/2(3X(mL - ml,,,,)- 1 « 9.2 to 121] 



arctan 



'IS 



(m. 



VNS' 



2m 2 K + m 2 )(m 2 



m: 



,1/2 



{2m 2 K — ml — m 2 )(m 2 , — m 2 ) 



41.9° 



(12) 



This kinematical QCD mixing angle ( |T2j ) or 9 = <fi — 54.74° w —12.8° has dynamical analogs 
p~6|] , namely the coupled SD-BS approach mentioned in the Introduction and used in Sec. 
IV below. Since this predicted 77-77' mixing angle in (|12D is compatible with the empirical 
values in (Rl) and (|9p, we use (fl~2| ) in the mixing angle relations (||) to infer the nonstrange 
and strange 77 masses, 



m ms = cos2 ^ m v + sin2 ^ m v' ~ ( 758 MeV )' 
= sin 2 <p m 2 + cos 2 m 2 , w (801 MeV) 2 



m 



(13a) 
(13b) 



Thus it is clear that the true physical masses 77(547) and 77'(958) are respectively much 
closer to the Nambu-Goldstone (NG) octet t/ 8 (567) and the non-NG singlet t/ (947) configu- 
rations than to the nonstrange 7/^(758) and strange 775(801) configurations inferred in Eqs. 
(|j~3|). However, the mean 77-77' mass (548 + 958)/2 w 753 MeV is quite near the nonstrange 
^ats(758). But since 77 8 (567) appears far from the NG massless limit we must ask: how close 
is ?7o(947) to the chiral-limiting nonvanishing singlet 77 mass? 

To answer this latter question, return to Fig. 1 and the quark annihilation strength 
(3 ~ 0.28 GeV 2 in Eq. (|TT|). These qq states presumably hadronize into the U^(l) singlet 
state I7/0) = \uu + dd+ss) /a/3, for effective squared mass in the SU(3) limit with (3 remaining 



unchanged ||34j| : 



m 



3(3 w (917 MeV) s 



(14) 
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This latter CL rjo mass in ([14]) is only 3% shy of the exact chiral-broken 770 (947) mass found 
in Eq. (|l]). (Such a 3% CL reduction also holds for the pion decay constant f n ~ 93 MeV 
-> 90 MeV H and for / + (0) = 1 -> 0.97 |§, the K-vr JQ 3 form factor.) 

Thus this 77-77' mixing resolution of the first Ua(1) problem is that the physical 7/(547) 
is 97% of the chiral-broken NG boson 7/3 (567). Also the mixing-induced CL singlet mass 
of 917 MeV in fllil) is 97% of the chiral-broken singlet 77 (947) in (H), which in turn is 99% 
of the physical 77' mass 7/ (958). This speaks to Weinberg's question |2j as to why there is 
no isoscalar, pseudoscalar Goldstone boson (with mass less than about \p3m^ ~ 240 MeV), 
associated with the spontaneous breakdown of the axial U^(l) symmetry. 



III. HADRONIC ETA DECAYS AND THE U A (1) PROBLEM 



As for the second Ua(1) problem, Weinberg in [Q] correctly identified the rapidly varying 
77 and 7T° poles for 77 — > 3tt° decay. However, one must also fold in the PCAC consistency 
approach of Refs. fll3l , |37| leading to the r\ — > 3ti amplitude magnitude with f n w 93 MeV, 



D7T 



Hera \v) 



(3/2#) tt u H em \t,) + OK/roJ) 



(15a) 



Here the factor of 3//^ on the RHS of ( 15a| ) corresponds to the three successive double pion 
PCAC reductions, while the factor of 1/2 characterizes Weinberg's || rapidly varying 77 and 



7r pole terms. Also this AI = 1 77 — > tt non-strong transition in (|15a|) reduces to [0,|T 



7T 



H em \V) 



COS< 



(tt° U 3 \l] NS ) 



cos42°(Am^ 



Ami 



-3900 MeV . 



(15b) 



In ( |15b|) we have invoked the CG u% = q\ 3 q quark tadpole (which is known []TT1 , |T2"1 to explain 
all P, V, B, D hadron SU(2) electromagnetic (em) mass splittings) using the SU(3) form 
(7r°|n 3 \rj NS ) = Am 2 K — Ami ~ —0.0052 GeV 2 , where Am 



K 



( 15b ) we have again invoked the 77-77' mixing relations 
(0)- 

Substituting Ql5b|) into (pa]) , one obtains the 77 — > 3n amplitude 



mj <+ — m^o, etc. Also in 



with mixing angle predicted by 



3vr l 



H em \ v ) = (3/2/£ 



7T 



0.68 



(16a) 



As for the experimental i] 3n o decay amplitude, taking a constant matrix element ( |16a| ) inte- 
grated over the Dalitz plot, one predicts an 77 — ► 3tt° decay rate 



T{r] 3n o) = (816 eV) 



3tt 1 



Hem \V) 



377 eV . (16b) 
rate of 380 ± 36 eV at the central 



The latter almost perfectly matches the 1998 PDG 
value. 

Alternatively we can extract the effective constant 3-body matrix elements A a , At,, A 
from data IHI 



r(rf 



3tt°) w 0.82 \AJ 2 keV 



► 37T U ) P 

tjttV) 



5.58 |AI keV 



0.38 keV, 
0.31 keV, 



1.06 \A r \ keV ss 42 keV, 



(17a) 
(17b) 
(17c) 
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leading to the dimensionless 3-body amplitudes 



\A a \ « 



0.68, \A h 



0.24, 



\AA « 



6.3. 



(17d) 



Note that the PCAC amplitude for (3tt°| H em \rj) in (pal) recovers the observed 77 — > 37r° rate 
in ( p. 6 b| ) or equivalently the constant Dalitz plot amplitude forms in ([17]) give |A | ~ 0.68 
which was earlier used to predict the 77 — > 37r° rate in Eqs. (|T6|) . 

This consistency pattern can also be applied to rj' — > 3n° decay presumably dominated 
by [38[] v( — > ?77r 7r followed by an em transition (7r°| i? em (77) : 



3tT 1 



#em W) 



71 



H em \ V ) (vrVr/IV) 



(m. 



3(3900 MeV 2 ) (6.3) (281000 MeV 



2\-l 



0.26. 



18a) 



In ( |18a| ) we have again used the em scale ( |15b|) (three times), the 77 propagator on the 7r° 
mass shell and the constant amplitude \A C \ ~ 6.3 in ( |1 7d| ) . The result 0.26 is near the 
constant amplitude \A b \ « 0.24 in ( |17d| ), or equivalently the i]' 3n decay rate is predicted to 
be 



r(r/ -> 3tt°) 



5.58 



3tt 1 



H em \rf) 



keV w 377 eV, 



;isb) 



near data § 313 ± 58 eV. 

Finally we consider the strong decays rf — > 7771 7r, with the charged to neutral pion branch- 
ing ratio being |J about 2, as expected via SU(2) symmetry. At first these decays were 
thought to be controlling the 77-77' mixing angle. Now, however, one begins by assuming an 
77-77' mixing angle [such as ~ 42° or 9 ~ —13° found earlier in Eqs. (|8|j9| JT2]) 1 , and then 
attempts to explain the observed 77' — > rj-n-ii rate given in Sec. I. 

To this end Singh and Pasupathy in Ref. |^] studied the 5 = a (983) scalar meson pole 
amplitude in 7/ — > 5tt, 5 — > 7/7r. Later Deshpande and Truong in also included a scalar 
meson a pole in this analysis with 77' — ► 770", cr — > 7r7r. These second authors in justified 
introducing this latter a in order to mask a soft-pion Adler zero which would drastically 
alter the tttt phase space. In fact the 77' — > 777r°7r° data shows only a small deviation from 
phase space, with linear amplitude A(l + ay) now requiring [§] a — —0.058 ± 0.013, and 
a = —0.08 ± 0.03 for 7/ — > i]7t + tt^ decay. 

Keeping only these two S and a pole terms, we slightly modify Refs. [3{J and write this 
combined 7/ — > riix + ix~ amplitude magnitude as 



.4 



^(77' — > 777r + 7T ) 



9Sr}w9ri' Sir 



+ 



ml — u — imsTs ' tti^ — s — imo-ro- 



(19) 



Here the combined 5 and a pole amplitudes have the same structure as in Ref. []39| except we 
always (rather than partially) keep the non-narrow widths || ~ 100 MeV and IV ~ 700 
MeV IRIIOfl- Also to estimate the pole denominators in flT5|), we follow Ref. [55] and take 



77?J — u ~ 2m v 'Ei w 2m r) i{m T) i — m,,) in the 77' rest frame with w w soft and 
s = [6.77-2.4t/] m2. 

Finally we choose the nonstrange a mass from the recent data analysis of Ref. [41] 



400 to 900 MeV 



mean mass m n w 650 MeV 



(20) 
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This is near e (700) used in [3J| and is supported by the 1998 PDG tables |J. Moreover a 
cr(650) is generated from linear a model (LcxM) dynamics [42] with LcrM coupling constants 
using the mixing relations ([]): 

/ m\ — vn? , \ 

9s m = cos <pgs VNS7r = cos0 I — INS I « 1.56 GeV, (21a) 

^ = sin 4>g Sms7r = sin j ^"^ 5 ) « 1.40 GeV, (21b) 

<W = ^ w /2/ w « 2.27 GeV, (21c) 
^jy'jjo- = cos (j) sin </><7a7T7r « 1.13 GeV . (21d) 

Note that gs VNS7T = g^mr in the chiral limit and also that the 77-77' mixing angle used (0 = 
41.9°) is as found from Eq. (0). 

Substituting the above numerical values back into ( |19"D leads to the 7/ — > r]7T + 7T~ ampli- 
tude magnitude 



2.20 2.57 
+ 



0.79-zO.lO 0.29-i0.46 



6.64 . (22) 



This LaM prediction in (|22| ) should be compared with the original estimates in |39| of 
\A\ m 8.5, a m —0.012. Also, \A\ m 6.64 in (|22]) is near \A C \ w 6.3 in (|17d|) assuming a 
constant matrix element and isospin invariance. Lastly accounting for the 7/ — > f]TT°TT° as 
well as the 77' — > 7]tx + tx~ amplitude, and folding in the slight Dalitz plot slope we predict the 
total decay rate (for the average slope a ~ —0.07): 

r(7/ -> 777T7T) = 3 \A\ 2 (1 + 0.24a + 0.27« 2 ) keV (23a) 
w 130 keV . (23b) 

This prediction ( |23b|) is in very good agreement with present data (131 ± 8 keV) as given 
in Sec. I. 

We differ from Ref. B^] primarily in that we use the LaM meson-meson couplings in Eqs. 



(pi|). An extraction of the 5t]ti coupling from the width of T(c'777r) ~ 100 MeV ^3[ gives for 
q = 321 MeV: 

T(5rrr) = ifejL , or \ 9&vr \ ~ 1.38 GeV . (24) 

87T7TT,| 

The latter coupling in (p4|) is reasonably near the LcxM coupling 1.56 GeV in (21a). 



IV. CONSISTENCY WITH DYNAMICAL CALCULATIONS 

As pointed out in the Introduction and Sec. II, there is a dynamical approach to the 



question of the Goldstone boson structure of the mixed 77(547) and t/(958) mesons [16 
namely the coupled SD-BS approach incorporating some crucial features of QCD, which 
leads to the similar conclusions on the mixing angle and masses as the analysis in Sec. II. 
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Before addressing its mass matrix, let us see what this approach tells us about the mixing 
angle that can be inferred from 77 decays. Since the SD-BS approach incorporates the 
correct chiral symmetry behavior thanks to D%SB and is consistent with current algebra, 
it reproduces (when care is taken to preserve the vector Ward-Takahashi identity of QED) 
the Abelian axial anomaly results, which are otherwise notoriously difficult to reproduce 



in bound-state approaches, as discussed in Ref. ||17|| . This gives particular weight to the 



constraints placed on the mixing angle 9 by the SD-BS results on 77 decays of pseudoscalars. 

A. 77 decays of the bound-state 7r°,r?, 7/ 

We express the broken-SU(3) pseudoscalar states vr°,?7 8 and r] through the quark basis 
states \ff) by 

\P) = E (^) l//>, (f = u,d,s), (25) 
where P = tt°, rj 8 , rjo simultaneously have the meaning of the respective indices j = 3, 8, 



on the SU(3) Gell-Mann matrices A J (j = 1,...,8) and on A = (y2/3)l 3 . This picks out 
the diagonal A 3 , A 8 , A in Eq. (|25|) . For future convenience we write the P(p) — > l{k)~f(k') 
amplitudes as 

T P (k 2 , k' 2 ) = £ ^ Q) f f f(k 2 , k' 2 ) , (26) 

where Tff(k 2 ,k' 2 ) = Tfj(k 2 ,k' 2 )/Q 2 are the "reduced" two-photon amplitudes obtained by 
removing the squared charge factors Q 2 from Tfj, the 77 amplitude of the pseudoscalar 
quark- ant iquark bound state of the hidden flavor //. 

The decay amplitudes (into real photons, k 2 = k' 2 = 0) of the physical states 77 and 7/, 
are given in terms of the predicted [16| 77 decay amplitudes of the SU(3) states r/ 8 and rj 
as 

r„(0,0) = cos 0^(0,0) -sin 0^(0,0) , (27) 
TV (0,0) = sin^ T m (0,0) + cos0r„ o (0,0) . (28) 



The best fit to the experimental 77 decay amplitudes was found in Ref. [16] for 9 — —12° 



for the concrete SD-BS model and parameters |L4| adopted there. In order to show that 
in the SD-BS approach 77 decays imply 9 somewhere in that ballpark (i.e., less negative 
than values favored by %PT) regardless of any model choice, and to be able to compare 
with other theoretical approaches which usually try to express P — > 77 amplitudes in terms 
of the leptonic (axial-current) decay constants fp, let us start with the light u, d sector in 
the chiral (and soft) limit. There, the SD-BS approach yields analytically and exactlyQ, and 
independently of the internal bound-state pion structure, 



1 The same holds [ 44 , 45 1 for the related process 7 — > tt + tt°tt . 
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7X0, 0) = 7^(0,0)= 7X0,0) 



N, 



(29) 



t;o(o,o) 



2^2^ y 




(30) 



Of course, the calculated fl^JlqJTql value of /„. does depend on the (modeling of the) internal 
pion structure, but the empirically successful axial-anomaly chiral-limit relation fl3"0D does 
not. 

The 7T° — > 77 decay amplitude for a possibly nonvanishing pion mass, can be used as a 
definition of pionic 77-decay constant f n by demanding that this amplitude be written in 
the form of the massless, CL amplitude (|30|), but with f n in place of /„.: 7X0, 0) = I/^tt 2 ^. 
Obviously, f n = f n in the CL, and f n is a convenient way to re-express the 77 amplitude 
in the case of a nonvanishing pion mass, because the Veltman-Sutherland theorem, PCAC, 
and the empirical success of the chiral-limit anomaly result fl3!]|), guarantee that f n ~ f n 
always holds for any realistic description of the light u, d sector. For simplicity of discussion, 
we therefore use f v = f n in this subsection, as the Veltman-Sutherland theorem guarantees 
that this can be wrong only by several percent. Although the chiral limit formula (|30| ) can 
be applied without reservations only to pions, it is customary to write the amplitudes for 
V8,Vo 77 hi the same form as (|3~0"D, defining thereby the 77-decay constants f m and f vo : 



7X0,0) 




U T^o(0,0) 



(31) 




fvo 



8TX0,0) 
73 ' 



(32) 



As pointed out by j46|, f m and f m 
current decay constants f m and /, 

through the 77-decay constants f m and f vo , yields the customary (see, e.g. [j46| ) forms for 
the rj and rj' decay widths: 



are not a priori simply connected with the usual axial- 
m , in contrast to f n ~ /„-. Expressing T m (0, 0) and T m (0, 0) 

e.g. 



r(r) 
T(rj' 



77) 
77) 



64vr 3 3/2 
64vr 3 3/2 



h 

fm 



cos 9 — V8 J ^r- sin 9 



Y~ sin 9 + 

JV8 



fvo 

7 

fvo 



fn ,1 

cos 9 



1 2 



(33) 
(34) 



The even more customary version of fl33|) and (p4|) in which the axial-current decay 
constants f V8 and f m appear in place of f m and f m requires a derivation where PCAC and 
soft meson technique are applied to the 77 — 77' complex [47j. For the indeed light pion, these 
assumptions are impeccable (leading to f n = f n ), but not for the 77-77' complex. For such 
a heavy particle as 77' they are quite dubious. However, we do not need and do not use 
these assumptions since we directly calculated the 773 and rj decay amplitudes, i.e., f m and 
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f vo , just as the axial-current pseudoscalar decay constants f V8 and f m were calculated fT6| 
independently of the 77 processes. In contrast to f n = f n , f V8 and f m cannot be equated, 
as the difference between them was found to be quite important [[T 



The precise values of f m and f m are model dependent, but f V8 < w f n holds in this 
approachF] generally, i.e., independently of chosen model details, as long as the s-quark mass 
is realistically heavier than the u, <i-quark masses. To see this, let us start by noting that 
fm < /t * s equivalent to T m (0, 0) > T w o(Q, 0)/a/3, and since we can re-write Eq. (^) for i] 8 

as 



r w (o,o) = 5!^ + ~[r (iJ (o,o)-r rf (o,o)] , (35) 

the inequality f V8 < f„ is in our approach simply the consequence of the fact that the 
("reduced") 77-amplitude of the ss-pseudoscalar bound state, T sS , is smaller than the cor- 
responding non-strange 77-amplitude T d( j (= T ua = T T o in the isosymmetric limit), for any 
realistic relationship between the non-strange and much larger strange quark masses. This is 
the reason why in this approach one cannot fit well the experimental 77, 7/ — > 77 widths with 
the mixing angle as negative as in chiral perturbation theory descriptions (9 ~ —20°), but 
rather with 9 ~ —12°. This is easily understood, for example, with the help of Fig. 1. of Ball 



et al. \ Wjj , where the values of f V8{0) / fn consistent with experiment are given as a function of 
the mixing angle 9. Their curve shows that values f m / fn < 1 permit accurate reproduction 
of 77,77' — > 77 widths only for values less negative than —15°. [It does not matter that 
they in fact plotted f Vs r J fn and not f V8 r J fn- Namely, they used Eqs. d33|)-(p4|) above for 
comparison with the experimental 77-widths, just with f m{0) / f n instead of f V8(0) / fn, so that 
the experimental constraints displayed in their Fig. 1 apply to whatever ratios are used 
in these expressions. One should also note that since in our approach f V8 , f Vo and are 
not free parameters but predicted quantities, the two widths 77, 77' — > 77 cannot be fitted 
exactly by adjusting just one parameter, 9. Rather, we fix 9 by performing a x 2 fit to the 
widths.] On the other hand, the more negative values 9 < —20° give good 77, 77' — > 77 widths 
in conjunction with the ratio f m l fn = f^lfn = 1-25 obtained by in x?T. However, 
the coupled SD-BS approach belongs among constituent quark approaches^ and for them, 
considerably less negative angles, 9 w —14° ± 2 [p9], are natural. 



Ref. |16| showed that these bounds and estimates are very robust under SD-BS model 
variations and can be taken as model independent. For example, for chiral u, d quarks, 



2 This is different from chiral perturbation theory fl"8fl . Nevertheless, for the axial-current decay 
constants our approach gives f V8 > f n (see Appendix or Ref. |]l6|). Even the numerical value 



obtained in our concrete SD-BS calculation [16|, f V8 = 1.31f n , is rather close to f V8 = 1.25/^ 
obtained in chiral perturbation theory 



3 However, at least one effective-meson-Lagrangian approach, that of Benayoun et al. [pQ] , yields 
results quite close to ours: 6 = —11.59° ± 0.76° and their Eq. (29), where their (/s, /1) correspond 
to our (f m J m ), with f 8 /U = 0.82 ± 0.02 and hi fn = 1-15 ± 0.02. 
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^ f ^ f 

fm = 5 - ' /%= 5 + ^f ss -(M)' (36) 

leading to the bounds | /„• < < /„- and /„- < f Vo < Also, considerations based 

on the Goldberger-Treiman relation showed that T S s(0,0) < T u u{0,0) is simply due to 
fss ~ fir+2(fK+ — fir) > fir (where f S g is the axial-current decay constant of the unphysical ss 
pseudoscalar bound state), and that a good estimate of the 77-amplitude ratio is the inverse 
ratio of the pertinent constituent quark masses: T sS (0, 0)/T nS (0, 0) ~ m/m s . Equations ( |36|) 
then give the relations [reducing to f m = and f no = f v in the U(3) limit, just like Eqs. 
(^) themselves] 

f ~ ^JjL f ~ ^JjL (-17) 

b — 2m/m s 5 + m/m s 



obtained also by Ref. |5T] using the simple quark loop model with constant constituent 
masses. These estimates are (for reasonable m/m s ) close to what Ref. []TB[ calculated with 
a concrete SD-BS model choice fll4j] , namely f Vo / f n = 1.067 and f m / f w = 0.797. For these 



concrete model values, 77,77' — > 77 widths (p3|)-(M) fit the data best for 8 = —12.0°. 



B. Introducing X into the SD-BS mass matrix 

For the very predictive SD-BS approach to be consistent, the above mixing angle ex- 
tracted from 77, 7/ — > 77 widths, should be close to the angle 9 predicted by diagonalizing 
the 77-7/' mass matrix. In this subsection, it is given in the quark // basis: 



M 2 = diag(M 2 „ M 2 ,-, Ml 



+ P 



1 1 1 
1 1 1 
1 1 1 



(38) 



As in Sec. II, 3/3 (called in Ref. ]IB[) is the contribution of the gluon axial anomaly to 
m 2 , the squared mass of 770. We denote by Mjj, the masses obtained as eigenvalues of the 
BS equations for qq pseudoscalars with the flavor content //' (/, /' = u,d,s). However, 
since Ref. |16| had to employ a rainbow-ladder approximation (albeit the improved one of 
Ref. Jy]]); it could not calculate the gluon axial anomaly contribution 3/3. It could only 
avoid the U,4(l)-problem in the 77-77' complex by parameterizing 3/3, namely that part of 
the 770 mass squared which remains nonvanishing in the CL. Because of the rainbow-ladder 
approximation (which does not contain even the simplest annihilation graph - Fig. 1), the qq 
pseudoscalar masses M^j, do not contain any contribution from 3/3, unlike the nonstrange 
and strange 77 masses m rjNS [in Eq. ( |13a| )1 and m Vs [in Eq. ( |13b| )1, which do, and which 
must not be confused with M uU = M d( j and M sS . Since the flavor singlet gluon anomaly 
contribution 3/3 does not influence the masses m n and rrix of the non-singlet pion and kaon, 
the realistic rainbow-ladder modeling aims directly at reproducing the empirical values of 
these masses: M uU = M d ^ = m v and M S( j = mj-. In contrast, the masses of the physical 
etas, m v and m v i, must be obtained by diagonalizing the r] 8 -r] sub- matrix containing both 
Mfj and the gluon anomaly contribution to m r 
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Since the gluon anomaly contribution 3(3 vanishes in the large N c limit as 1/N C , while 
all Mfj, vanish in CL, our qq bound-state pseudoscalar mesons behave in the N c — > oo and 
chiral limits in agreement with QCD and xPT (e.g., see []52|): as the strict CL is approached 
for all three flavors, the SU(3) octet pseudoscalars including rj become massless Goldstone 
bosons, whereas the chiral-limit-nonvanishing //-mass 3(3 is of order 1/N C since it is purely 
due to the gluon anomaly. If one lets 3(3 — ► (as the gluon anomaly contribution behaves 
for N c — ► oo), then for any quark masses and resulting Mtt masses, the "ideal" mixing 
{6 = —54.74°) takes place so that 77 consists of u, d quarks only and becomes degenerate 
with Ti, whereas 77' is the pure ss pseudoscalar bound state with the mass M S g. 

In Ref. ||16|| , numerical calculations of the mass matrix were performed for the realistic 
chiral and SU(3) symmetry breaking, with the finite quark masses (and thus also the finite 
BS qq bound-state pseudoscalar masses Mfj) fixed by the fit [14j] to static properties of many 
mesons but excluding the 77-77' complex. The mixing angle which diagonalizes the 773-770 mass 
matrix thus depended in Ref. [IB] only on the value of the additionally introduced "gluon 
anomaly parameter" 3(3. Its preferred value turned out to be 3(3 = 1.165 GeV 2 =(1079 
MeV) 2 , leading to the mixing angle 6 = —12.7° [compatible with (j) = 41.9° in Eq. (|i~2|)1 
and acceptable rj — > 77 and 77' — > 77 decay amplitudes. Also, the 77 mass was then fitted 
to its experimental value, but such a high value of 3(3 inevitably resulted in a too high 77' 
mass, above 1 GeV. (Conversely, lowering 3(3 aimed to reduce m v >, would push 9 close to 
—20°, making predictions for 77,77' — > 77 intolerably bad.) However, unlike Eq. flip] ) in the 
present paper, it should be noted that Ref. |16j did not introduce into the mass matrix the 
"strangeness attenuation parameter" X which should suppress the nonperturbative quark 
// — > /'/' annihilation amplitude (illustrated by the "diamond" graph in Fig. 1) when / 
or /' are strange. 

On the other hand, the influence of this suppression should be substantial, since X m 
rh/m s should be a reasonable estimate of it, and this nonstrange-to-strange constituent mass 
ratio in the considered variant of the SD-BS approach |T5| is not far from X in Eq. (|TI|) and 
from the mass ratios in Refs. |31|-p3|, and is even closer to the mass ratios in the Refs. ||30|| . 
Namely, two of us found |1| it to be around M u (0)/ Ai s (0) = 0.615 if the constituent mass 
was defined at the vanishing argument q 2 of the momentum-dependent SD mass function 

mm 2 )- 

We therefore introduce the suppression parameter X the same way as in the NS-S mass 
matrix ([TOD, whereby the mass matrix in the // basis becomes 



M 2 = dmg(M 2 a ,M 2 (I ,M 2 ,) 



1 1 X 
1 1 X 

XXX 2 



(39) 



In a very good approximation, Eq. ( |39| ) recovers (in the tt°-NS-S basis) Eq. ( |T0| ) for the 
2x2 77-77' subspace. This is because M 2 S differs from 277^ — m 2 only by a couple of percent, 
thanks to the good chiral behavior of the masses M^p calculated in SD-BS approach. (These 
M 2 j, and the CL model values of f n and quark condensate, satisfy Gell-Mann-Oakes-Renner 
relation to the first order in the explicit chiral symmetry breaking ||15|| .) The SD-BS- 
predicted octet (quasi-) Goldstone masses Mjj, are known to be empirically successful in our 
concrete model choice fT4]| , but the question is whether the SD-BS approach can also give 
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some information on the X-parameter. If we treat both 3(3 and X as free parameters, we 
can of course fit both the rj mass and the rj mass to their experimental values. For the 
model parameters as in Ref. |14| (for these parameters our independent calculation gives 
= M uU = 140.4 MeV and M sS = 721.4 MeV), this happens at 3(3 = 0.753 GeV 2 =(868 
MeV) 2 and X = 0.835. However, the mixing angle then comes out as 9 = —17.9°, which is 
too negative to allow consistency of the empirically found two-photon decay amplitudes of 
r) and rj', with predictions of our SD-BS approach for the two-photon decay amplitudes of 
i] 8 and rjo [16] . 



Therefore, and also to avoid introducing another free parameter in addition to 3(3, we 
take the path where the dynamical information from our SD-BS approach is used to estimate 
X. Namely, our 77 decay amplitudes Tfj can be taken as a serious guide for estimating the 
X-parameter instead of allowing it to be free. We did point out in Sec. II that the attempted 
treatment |34| of the gluon anomaly contribution through just the "diamond diagram" con- 



tribution to 3(3, indicated that just this partial contribution is quite insufficient. This limits 
us to keeping 3(3 as a free parameter, but we can still suppose that this diagram can help 
us get the prediction of the strange-nonstrange ratio of the complete pertinent amplitudes 
// — ► /'/' as follows. Our SD-BS modeling in Ref. |16| employs an infrared-enhanced 
gluon propagator JT4] , |T7| ] weighting the integrand strongly for low gluon momenta squared. 
Therefore, in analogy with Eq. (4.12) of Kogut and Susskind || (see also Refs. |53| , |54]| ) , we 



can approximate the Fig. 1 amplitudes // — > 2gluons — > /'/', i.e., the contribution of the 
quark- gluon diamond graph to the element //' of the 3x3 mass matrix, by the factorized 
form 

r // (0,0)CT / , / ,(0,0). (40) 

In Eq. (fffi|), the quantity C is given by the integral over two gluon propagators remaining 
after factoring out 7^(0,0) and Ty,j,(0,0), the respective amplitudes for the transition of 
the qq pseudoscalar bound state for the quark flavor / and /' into two vector bosons, in this 
case into two gluons. The contribution of Fig. 1 is thereby expressed with the help of the 
(reduced) amplitudes 7//(0, 0) we calculated for the transition of qq pseudoscalars to two 
real photons (k 2 = k' 2 = 0). Although C is in principle computable, all this unfortunately 
does not amount to determining (3, (3X and (3X 2 in Eq. (|39| ) since the higher (four-gluon, 
six-gluon, ... , etc.) contributions are clearly lacking. We therefore must keep the total 
(light-)quark annihilation strength (3 as a free parameter. However, if we assume that the 
suppression of the diagrams with the strange quark in a loop is similar for all of them, 
Eq. (|40| ) and the "diamond" diagram in Fig. 1 help us to at least estimate the parameter X 
as X « T S s(0, 0)/T U u(0, 0). This is a natural way to build in the effects of the SU(3) flavor 
symmetry breaking in the qq annihilation graphs. 

We get X = 0.663 from the two-photon amplitudes we obtained in the chosen SD-BS 



model ||14|| . This value of X agrees well with the other way of estimating X, namely the 
nonstrange-to-strange constituent mass ratio of Refs. |]3~T|-|3"B|] . With X = 0.663, requiring 
that the 2x2 matrix trace, m 2 + m 2 ,, be fitted to its empirical value, fixes the chiral- 
limiting nonvanishing singlet mass squared to 3(3 = 0.832 GeV 2 =(912 MeV) 2 , just 0.5% 
below Eq. fll4|). The resulting mixing angle and rj, rf masses are 

9 = -13.4° , m v = 588 MeV , rry = 933 MeV . (41) 
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The above results of the SD-BS approach |RJ are very satisfactory since they agree well 
with what was found in Sec. II by different methods. Let us close this section by exploring 
the stability of these results on model variations. Except the introduction of 3/3(= A^), these 
SD-BS results were obtained without any other parameter fitting, with the model parameters 



resulting from the very broad previous fit [|14]], but actually giving us, in our independent 
calculation, a few percent too high results for m n and w,k- To possibly improve, and in 
any case check the robustness of the consistency with Sec. II (and subsection IV.A) on 
variations of our model description, we therefore perform a refitting in the sector of u, d and 
s quarks, to reproduce exactly the average isotriplet pion mass = M uU = 137.3 MeV 
and isodoublet kaon mass = 495.7 MeV. As Table I shows, the changes are small, and 
lead to Mu(0)/M s (0) = 0.622 and X = f sS (0, 0)/f uU (0, 0) = 0.673. Using this X to fit the 
sum of the squared rj and rf masses to the empirical value, yields the column B in Table II, 
where we see a slight improvement in the 77 and rf masses with respect to the results (fXTI), 
while the mixing angle is still acceptable, being less than 2° away from the angle favored in 
Sec. II. 

If we treat X as the second free parameter (this procedure yields the column C of Table 
II) so that we are able to fit and precisely to their experimental values, we get 
X = 0.805, along with the mixing angle 9 = —14.9° and the chiral-limit-nonvanishing 
singlet mass 3/3 = 0.801 GeV 2 =(895 MeV) 2 . This is noticeably closer to 9 and 3f3 resulting 
from other procedures (where X is not a free parameter) than before the aforementioned 
7T° - K refitting to m n = 137.3 MeV and m K = 495.7 MeV. 

Next, we note in the column D of Table II that the slightly improved fit to the masses 
also led to somewhat improved 77,77' — > 77 widths when we extract from them 9 = —12.8°, 



practically the same as Ref. |]T6| and the Sec. II result (y^)- All the three possibilities B, C, 
and D, do not differ too much from each other, and agree reasonably with the experimental 
masses and 77 widths given in column E as well as with the corresponding results of Sec. 
II. This contrasts with column A, which also contains the results of the new fit but with 
X = 1. Column A shows that when X = 1, a good description of the masses requires a 9 
value too negative for a good description of the 77 widths in the SD-BS approach. Column 
A thus convinces us that it was precisely the lack of the strangeness attenuation factor X 
that prevented Ref. |16| from satisfactorily reproducing the 77' mass when it successfully did 
so with the 77 mass and 77 widths. 

V. CONCLUSION 



In Sec. II we studied the first U^(l) problem associated with the Goldstone structure 
of r/(547) and 7/(958) mesons. Following a QCD gluon-mediated approach to 77-77' particle 
mixing, we began by extracting an 77-77' mixing angle <fi ~ 42° in the NS-S basis or 9 ~ —13° 
in the singlet-octet basis. This led to eta masses 773 (567), 770 (947) with chiral- limiting (CL) 
77o(917). Then the physical eta mass 7/(547) is 97% of 77 8 (567), while 7/(958) is 104% of 
the CL 77o(917). Such a 3-4% CL suppression is likewise found for the pion decay constant 
f n « 93 MeV 90 MeV and for the K l3 form factor / + (0) = 1 -> 0.96-0.97. 

Then in Sec. Ill we studied the second U^(l) problem associated with eta meson hadronic 
decay rates. The 77,77' — ► 3ti° (AJ = 1) decay rates of 377 eV followed from PCAC Con- 
sistency. Also a (strong) decay rate of 130 keV for 77' — > rjnn was obtained from S and a 
scalar meson poles combined with linear a model couplings. These three rates are com- 
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patible with data finding @ T(t] -> 3tt°) = 380 ± 36 eV, T(rj' -> 3tt°) = 313 ± 58 eV and 
r(r/ -> 777T7r) = 131 ± 8 keV. 

Finally, in Sec. IV we showed the consistency of the above results with those obtained in 
a chirally well-behaved quark model which was explicitly constructed through D^SB, SD and 
BS equations. For example, described variations of our SD-BS approach lead to 9 ~ — 13°±2° 
and to the corresponding CL rjo mass \/3/3 = 912 ± 18 MeV. Successful reproduction of the 
Abelian axial anomaly amplitudes in the CL in this bound-state approach, gives particular 
weight to our conclusion that so far away from the CL as in the case of the rj-rf complex, 
77-decay constants (/^ g , f Vo ) differ significantly from the usual axial-current decay constants 
(fvsifrio)- By allowing for the effects of the SU(3) flavor symmetry breaking also in qq 
annihilation graphs, we have improved the 77-?/ mass matrix with respect to the mass matrix 
in Ref. |16| (via the strangeness attenuation factor X = 0.663). 



The consistency of our approach and results with the two-mixing-angle scheme is shown 
in detail in the Appendix, where we also compute the mixing angles and axial decay constants 
in that scheme. 
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APPENDIX: CONNECTION WITH THE TWO MIXING ANGLES SCHEME 

In this appendix, we clarify the relationship of our approach with the two-mixing-angle 
scheme considered by Leutwyler and Kaiser ]2"2"| , P3"| as well as FKS , and reviewed by 

Feldmann |28| . 



The two-mixing-angle scheme is defined with respect to the mixing of the decay constants. 
As such, it is very suitable in studies where manipulating decay constants is crucial, e.g., 
when one expresses amplitudes through them with the help of PCAC. In other situations, 
the mixing of the states may be crucial. The SD-BS variant of our approach, where one 
explicitly solves for quark-antiquark bound states, and then uses them for direct calculation 
of amplitudes, is an especially clear example of that. (See Ref. [[^| for another recent 
example.) As FKS themselves state several lines below their Eq. (1.3), the appearance 



of the four parameters in the two-mixing-angle scheme, namely /§, fo, 9$ and 60, raises anew 
the problem of their mutual relations and their connection with the mixing angle of the 
particle states, which is necessarily a single one. In our approach, it is convenient to utilize a 
mixing scheme defined with respect to a state basis corresponding to the broken SU(3) flavor 
symmetry, \r] NS ) and \i] S ) or, equivalently, the effective?] SU(3)-broken \rj 8 ) and \r] ) [Eqs. 
(H)], i.e., the state mixing angle of Eqs. (B) or (mathematically completely equivalently) 



4 Note that in spite of the differences in notation, the effective octet and singlet states (||) of the 
broken SU(3) flavor, correspond to the effective octet and singlet states 77s and 770 given, e.g., by 
Eq. (85) in Feldmann's review |28||. 
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the state mixing angle 9 of Eqs. (|2|). Nevertheless, we show below that a) we can calculate 
quantities utilized in the two-mixing-angle scheme defined with respect to the mixing of the 
decay constants, and b) what we find for these quantities is close to what is quoted in Refs. 



Independently of any specific approach, one can always define quite generally the axial- 
current decay constants /®, f*,, /°, and /£ as the matrix elements 



(Q\A a »(x)\Ptp))=ift 1 fe 



ip-x 



0; P 



(Al) 



These definitions are to be contrasted with the somewhat arbitrary definitions of the just 
two individual axial-current decay constants f v and f v r often used in this context, where 
f v = cos 2 /ivs+sin 2 fs and = sin 2 <fi Jns+cos 2 <fi fs, where f^s an d fs are defined below 



in Eqs. ( |A6| ) . These f v and f v > stem from Eqs. (]£]) in conjunction with the rather arbitrary 
definitions AHx) = cos0 A^ s (x) — sin0 Ag(x) and AZ(x) = sin0 A^ s (x) + cos0 A$(x), 
with Eqs. (|A3|) below defining the nonstrange and strange axial currents, A^ s (x) and Ag(x). 
In contrast to this, the four decay constants /®, /°, and /°/ are precisely defined by Eqs. 
( |AT] ) and therefore can have unambiguous and process- independent meaning ||28|| . 

Following the convention of Leutwyler and Kaiser [^,^], the four decay constants 
frfi /°; an d /n' are parametrized in terms of two decay constants f , / 8 , and two angles 9 , 
8 : 



We define the currents: 
1 



A% s (x) 



A Sli (x) 



-A^\ 



x 



f 

Jrf 

f 

fO 
Jrf 



'-A Q n 



cos 9 8 f 8 , 
sin 6» 8 f 8 , 
- sin 6» f 
cos 6»o f ■ 



u(x)^^u(x) + (i(x)7 M 75(i(a; 



V3 



A 0fl ( 



V2 



s(x)7 M 7 5 s(x) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 



(A3a) 
(A3b) 



The corresponding NS-S decay constants (analogous to the constants /5, /°, and /S 
defined above) are defined as 



(0|^(x)|P(p))=^p^e 



-ip-x 



F = NS,S; P = r),r)' 



(A4) 



The relations ( |A3| ) between the currents dictate that the relations between these two sets of 
decay constants are given, exactly and model independently, by 



fNS fS 

J ri J r> 

fNS fS 

Jrf Jrf 



f8 fO 

J r\ J rj 

f8 fO 

Jrf Jrf 



1 

3 



1 



(A5) 



where we used matrix notation for compactness. 
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If we have well-defined nonstrange— strange states \r] N g) and \rjs) (§) [as in our SD-BS 
approach] we can define the decay constants f^s and fg through 



(0\A% s (x)\r ]N s(p))=if N sp^ i ^ 
(0\A» s (x)\r )s (p))=ifsP tl e- i *>-* , 
(0\A% s (x)\ Vs (p)) = , 
(0\A^x)\ VNS (p)) = . 



(A6a) 
(A6b) 
(A6c) 
(A6d) 



Since the states \rj) and \rf) are given by Eqs. (f|) as the linear combinations of \t]ns) and 
1 77s), we can relate the constants {f* s , f$ s , f^ } f%} with {f NS , f s }: 



fNS fS 

J ri J n 

fNS fS 

Jn' Jn 1 



COS( 

sin ( 



• sin i 

COS) 



fNS 

fs 



(A7) 



Using Eq. 



we can relate the decay constants {/®, f*,, /°, /£} with {fNS, fs} 



fS f0~ 

J r) J r\ 




f8 fO 
_ J v ' Jrf _ 





cos (p — sin i 
sin 6 cos i 



fNS 

fs 



1 2 
V 3 
F2 J_ 

3 y/3 



(A8) 



which is the samef] as in Feldmann et al. |[26|| , who use the notation / g = f^g, f s = fg. 

The above equations together with the definitions ( |A2|) give the following solutions for 
fs, fa, 9s, and O : 




1 2 



arctan f #1 

V JNS J 



2 1 



— arctan 



V2f 



NS 



fs 



(A9a) 
(A9b) 

(A9c) 
(A9d) 



These relations (the same as already found by FKS [2J-g8|) are obtained in a general way, 
independently of specifics of any given approach; therefore, we can also apply them within 
our framework. In Sec. II, we have already pointed out the agreement of our results for 
our preferred state mixing angle <fi ~ 42° with the FKS results quoted in Refs. [24-28], 
but the values we find for /q, /§, 9q, and 8s are also similar to theirs. The important 
quantity here is y — fNs/fs, giving the extent of the SU(3) breaking. In our approach, 
this quantity is also present and plays a crucial role. The information it carries enables 



5 Note that the last matrix in our Eq. ( |A8| ) is just W {Oideal) i n the notation of Ref. [26], where 
Oideal = arctan V2. 
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us to calculate what the decay-constant-mixing angles 9g and 9q would be in our approach 
(although we stress again that because of the way we calculate, their "rough average" 8 
must retain the central role because it has the meaning of the state-mixing angle). This is 
all because the y-ratio is essentially (as easily seen through the Goldberger-Treiman relation 
for constituent quarks) our parameter X, which can be (at least approximately) expressed 
as the ratio of the nonstrange-to-strange constituent quark mass: X w rh/m s . In addition, 
in our very predictive coupled SD-BS approach, we can directly calculate all decay constants 
including f g = f N $ and f s = fs, and this again gives (in a good approximation) the same 
value for X — y — fNs/fs- While we have f^s — fn (also assumed by FKS in 
their theoretical analysis), our calculation yields fs = 1.4505/^, less than 3% more than 
the theoretical FKS prediction |26|-p8|. (Note that we used the symbol f sS for f s = fs 
Our chosen model therefore gives y = fNs/fs = 0.6894, leading to fs = 1.318/^ and 



fo = 1.170 f n . (Interestingly, this is practically equal to our SD-BS model values f m = 1.31/V 
and f vo = 1.16 f n [HJ for the octet and singlet axial-current decay constants f m and f Vo , 
mentioned in Sec. IV. They are defined in the standard way through the matrix elements 
(0\A a ^\r] a ) , (a = 8,0), so that the definitions (£|) and ( |A3| ) imply that f m and f m are 
straightforwardly expressed through f^s and fs by the relations f m = ^fNS + |/s an d 
fm = ^fxs+ ^fs- We thus note that the quadratic relations ( |A9a| ) and ( |A9c| ) for differently 
defined octet and singlet constants / 8 and f , lead to similar values as the linear relations 
for f m and f m .) 

Using in Eqs. ( |A9|) our preferred state mixing angle <fi = 42°, our model value y = 
fNs/fs = 0.6896 also leads to the following decay-constant-mixing angles in the r]$ — r] 
basis: 8$ = —22° and 9 = —2.3°, close to the theoretical FKS results [EME7|. See also 



Table 1 in Ref. line "FKS scheme & theory", giving the 9 8 = -21.0° and 9 = -2.7°, 
while the line "FKS scheme & phenomenology" in the same table has only somewhat more 
negative #o but larger fo/f-K- The "FKS scheme & theory" then implies the state-mixing 
angle 9 ~ —12°, in agreement with our results. This is as it should be, as we note that the 



mass matrix in Ref. p8 |, its Eq. (73), coincides with ours when the anomaly contribution 
a 2 is identified with the "Veneziano term" A^/3(= (3), as we do. Recalling that we have 
already pointed out the agreement of the NS-S mixing angles obtained by us and by FKS 
|28| , one can see that everything tallies. 
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TABLES 



p 


mp 


fp 


T P (0,0) 


M q (0) 


IT 


0.1373 


0.0931 


0.257 


0.374 


K 


0.4957 


0.113 






ss 


0.7007 


0.135 


0.0815 


0.601 



TABLE I. The first column displays results of refitting tt, K and ss masses (m sS = M sS ) 
obtained in the qq bound state SD-BS approach with the slightly changed explicit chiral symmetry 
breaking bare masses fh u ^ = 2.965 MeV and rfi s = 69.25 MeV. These it and ss-pseudoscalar masses 
are input parameters for 77-77' fit in Table II. The last column is the constituent quark mass M q (0) 
pertinent to the corresponding qq meson, namely M u (0) = Md(0) for the pion and M. s (0) for the 
unphysical ss pseudoscalar. The masses mp and Ai q (0) as well as the pseudoscalar axial-current 
decay constants fp are in units of GeV, while the 77 decay amplitudes Tp(0, 0) are in GeV -1 . 





A 


B 


C 


D 


E 


X 


1.0 


0.673 


0.805 






3/3 


0.707 


0.865 


0.801 






e 


-19.5° 


-11.1° 


-14.9° 


-12.8° 




m v 


0.5048 


0.5777 


exp. 




0.54730 




0.9809 


0.9398 


exp. 




0.95778 


L(?7 77) 


0.63 


0.44 


0.52 


0.48 


0.46 ± 0.04 


IV - 77) 


3.61 


4.61 


4.16 


4.41 


4.26 ±0.19 



TABLE II. Various fits for the masses, mixing angle 9, and 77 decay widths L in the 77-77' 
complex. In columns A, B, and C, the free parameter in the mass matrix is 0, whereas X is free 
only in column C. Column A: results with X = 1 fixed by hand and fixed by fitting + m'^ 2 . 
Column B: results with X estimated from the ratio of the reduced strange and nonstrange 77 
amplitudes, and (3 fixed by fitting 777,^ + ttt,^ 2 . This column gives the best predictions, especially 
considering its only free parameter is 0. Column C: results with X treated as the second free 
parameter, making possible that 777,^ and 777,^/ are both fitted to their experimental values exactly. 
Column D: fitting the empirical 77 widths of 77 and 77' with 9 as the free parameter (and empirical 
777,7 and m v >), independently of the masses and 9 obtained from the mass matrix considerations. 
Column E: experimental values. Among the dimensionful quantities, 30 is in units of GeV 2 , m v 
and m' in GeV, while T(rj — > 77) and L(?7' — > 77) are in units of keV. 
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APPENDIX: FIGURE CAPTIONS 

Fig. 1: Nonperturbative QCD quark annihilation illustrated by the diagram with two-gluon 
exchange. It shows the transition of the // pseudoscalar P into the pseudoscalar P' 
having the flavor content /'/'. The dashed lines and full circles depict the qq bound- 
state pseudoscalars and vertices, respectively. 
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